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Abstract 

We analyze, in the paradigm of open quantum systems, the reduced dynamics of a freely falling 
two-level detector in de Sitter space-time in weak interaction with a reservoir of fluctuating quan- 
tized conformal scalar fields in the de Sitter-invariant vacuum. We find that the detector is asymp- 
totically driven to a thermal state at the Gibbons-Hawking temperature, regardless of its initial 
state. Our discussion therefore shows that the Gibbons-Hawking effect of de Sitter space-time 
can be understood as a manifestation of thermalization phenomena that involves decoherence and 
dissipation in open quantum systems. 
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The theory of quantum open systems has, in recent years, greatly advanced our under- 
standing of fundamental issues at the foundation of quantum mechanics and nonequilibrium 
statistical physics, and it has been fruitfully applied to the nascent yet fast-growing fields 
of quantum information science, modern quantum optics, atomic and many-body systems, 
soft condensed matter physics, and biophysics. In this Letter, we will apply the open quan- 
tum system approach to the investigation of the Gibbons- Hawking effect [l| in the hope of 
gaining insights into de Sitter space-time from the prospective of quantum open systems. 
Our interest in this issue also lies in the fact that the AdS/CFT correspondence between 
quantum gravity on an AdS space-time and a conformal quantum field theory without grav- 
ity on a boundary has provided new and rewarding lines of research into many branches of 
physics, such as condensed matter physics (see Ref. j^] for a recent review) and quantum 
chromodynamics (see Ref. for a recent review), and there may also exist a holographic 
duality between quantum gravity on de Sitter space-time and a conformal field theory living 
on the boundary identified with the timelike infinity of de Sitter space-time 0|. 

Our plan is to examine the time evolution of a freely falling two-level detector in interac- 
tion with fluctuating vacuum conformal scalar fields in de Sitter space-time. The detector is 
treated as an open quantum system and the vacuum with fluctuations of the quantum fields 
as the environment. The evolution of the detector is subject to the effects of decoherence 
and dissipation due to its interaction with the environment, and as for any open system, 
the full dynamics of the detector can be obtained from the complete time evolution describ- 
ing the total system (detector plus the environment) by integrating over the field degrees 
of freedom, which are in fact not observed. We will show that the detector in de Sitter 
space-time, regardless of its initial state, is asymptotically driven to a thermal state at the 
Gibbons- Hawking temperature. The open quantum system approach in our Letter there- 
fore demonstrates that the Gibbons-Hawking effect can be understood as a manifestation of 
thermalization phenomena in the framework of open quantum systems. It is worth noting 
here that an examination of similar issues, i.e., the Unruh effect and the Hawking radiation, 
in theparadigm of open quantum systems has been carried out, respectively, in Ref. Q and 
Ref. 

When vacuum fluctuations are concerned in a curved space-time, one first has to specify 
the vacuum states. The vacuum states in de Sitter space-time can be classified into two 
categories: one is the de Sitter-invariant states, the others are those which break de Sitter 
inwiance 0. Generally, the de Sitter-mvariant vacuum, whose status in de Sitter space- 
time is just like Minkowski vacuum in the flat space-time, is deemed to be a natural vacuum. 
So we will consider the evolution in the proper time of a freely falling detector in interaction 
with a quantized conformally coupled massless scalar field in the de Sitter-invariant vacuum. 

We assume the combined system (detector + external fluctuating vacuum fields) to be 
initially prepared in a factorized state, with the detector having two internal energy levels 
and the fields in the de Sitter vacuum. Thus the detector can be fully described in terms 
of a two-dimensional Hilbert space, so that its states can be represented by a 2 x 2 density 
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matrix p, which is Hermitian = p, and normahzed Tr(p) = 1 with det(p) > . Without 
loss of generahty, we take the total Hamiltonian for the complete system to have the form 

H = Hs + H^ + \H' . (1) 

Here Hg is the Hamiltonian of the detector, which is taken, for simplicity, to be 

H. = if , (2) 

where is the Pauli matrix and Uq the energy level spacing. is the standard Hamiltonian 
of conformal scalar fields in de Sitter space-time, details of which need not be specified here 
and H' is the interaction Hamiltonian of the detector with the external scalar fields and is 
assumed to be given by 

H' = a^^{x) . (3) 

Let us note here that we can also write Hg and H' in more general forms j^, but that 
does not change the result of this Letter. In order to achieve a rigorous, mathematically 
sound derivation of the reduced dynamics of the detector, we will assume that the interaction 
between the detector and the scalar fields is weak so that the finite-time evolution describing 
the dynamics of the detector takes the form of a one-parameter semigroup of completely 
positive maps jsl, [oj. It should be pointed out that the coupling constant A in ([T]) should 
be small, and this is required by our assumption that the interaction of the atom with the 
scalar fields is weak. 

Initially, the complete system is described by the total density matrix ptot = p(0) (g) |0) (0| , 
where p(0) is the initial reduced density matrix of the detector and |0) is the de Sitter- 
invariant vacuum state for field In the frame of the detector, the evolution in the 
proper time r of the total density matrix ptot of the complete system satisfies 

= -iLH[ptot{r)] , (4) 

where the symbol Lh represents the Liouville operator associated with H 

Lh[S]^[H,S]. (5) 

The dynamics of the detector can be obtained by tracing over the field degrees of freedom, 
i.e., by applying the trace projection to the total density matrix p(r) = Tr$[ptot(r)] . 

In the limit of weak coupling which we assume in this Letter, the reduced density is found 
to obey an equation in the Kossakowski-Lindblad form 10|, III 

M!:^ = -.[i/eff,p(r)]+/:[p(r)], (6) 

where 

1 ^ 

>C[p] = - ^ aij [2 (jjp ffi - (Ji(jj p- p o-jO-j] . (7) 
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The matrix aij and the effective Hamiltonian Hcs are determined by the Fourier trans- 
form ^(A) and Hilbert transform /C(A) of the field vacuum correlation functions (Wightman 
functions) 

G+(x-y) = (0|$(x)<l>(y)|0), (8) 

which are defined as 



dre'^^ G+{x{t)) 



/C(A) = — / du 



g(^) 



Then the coefficients of the Kossakowski matrix be written explicitly as 



(9) 
(10) 

(11) 



with 



A = ^[g{uo) + g{-coo)] , B = ]^[g{ujo) - g{-ujo)] , c = g{Q) - a . (12) 

The effective Hamiltonian i/efr contains a correction term, the so-called Lamb shift, and one 
can show that it can be obtained by replacing ujq in Hg with a renormalized energy level 
spacing VL as follows 



eff 



0-3 = {^0 + «[/C(-wo) - /C(wo)]}o-3 , 



(13) 



where a suitable subtraction is assumed in the definition of /C(— wq) — /C((X'o) to remove the 
logarithmic divergence which would otherwise be present. 

In order to find out how the reduced density evolves with proper time from Eq. dH]), 
we need the Wightman function for the conformally coupled scalar fields in the de Sitter- 
invariant vacuum. Let us note that there are several different coordinate systems that 
can be chosen to parametrize de Sitter space-time 12(]. Here we choose to work with the 
global coordinate system {t, x, d, 4>) in which the freely falling detector is comoving with the 
expansion. The line element is 



ds"^ = df - cosh^()f:/a)[c?x^ + sinx^(rf6'2 + sin^ Bd^p^)] 



(14) 



with a = 3^/^A~^/^, where A is the cosmological constant. The parameter t is often called the 
world or cosmic time. The canonical quantization of a massive scalar field with this metric 
has been done in Ref. 

0,[l2|-0- In coordinates dTl), the wave equation for a massive scalar 

field becomes 



1 


"1 


cosh^ - ^ ) 


cosh^ t/a 




V a dt J 



A 



a^cosh t/a 



(15) 



where the Laplacian 



A 



sin^X 



dx\^ ^dx) sin 



d 
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sin 9 dip"^ 



(16) 
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From ( !T5|) one can get the eigenmodes, and define a de Sitter-invariant vacuum. Then the 



Wightman function can be written as [l8| 



1 i-^' 



G+(x(r), x(r'))) = ^ F { I + u,l - 



IGvra^ cosTTz/ 

where F is a hypergeometric function and 

t t' t t' 

Z{x, x') = sinh — sinh cosh — cosh — cos Q 

a a a a 

cos Q = cos X cos x' + sin X sin x'[cos ^ cos 0' + sin 6 sin 6*' cos{ip — ip')]' 



(17) 



1/2 



(18) 



In the massless, conformal couphng hmit, for a freely falhng detector, the Wightman function 
becomes 

1 



G+(x(r),x(r'))) = - 
Its Fourier transform can be found to be 



IGrr^a^ sinh^(^^2^ — ie) 



X e 



27roA 



271 e 



2naX-l 



(19) 



(20) 



This leads to 



A = i/2[g,si\) + g,,i-\)] 



A coth(7raA) 
Air 



B = l/2[g,,i\)-g,s{-\)] = ^. (21) 

4:71 



In order to solve Eq. ([6]) to find out how the reduced density evolves with proper time, let 
us express it in terms of the Pauli matrices, 
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Pir) 



(22) 



Substituting Eq. fl22|) into Eq. (|6]), one can show that the Bloch vector |p(t)) of components 
{Pi{t),P2{t),P3{t)} obeys 



|:|p(r)) = -2H|p(r)) + |r^) 



(23) 



where {r]) denotes a constant vector {0, 0, —45}. The exact form of the matrix H reads 

^ 2A + C Q/2 ^ 
-n/2 2A + C 
2A 



\ 



(24) 



This matrix is nonsingular and its three eigenvalues are Ai = 2 A, A± = {2 A + C) ± iQ/2. 
Since the real parts of these eigenvalues are positive, at later times, |p(t)) will reach an 
equilibrium state \p{oo)) [2l|, which can be found by formally solving Eq. fl23|) 



\p{r)) 



'2Ht 



|P(0)) + (1 



-2Ht 



)\Poo) 



(25) 
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with 





VJ 



(26) 



If we let P = 1/T = 2 aicta.nh.{B /A)/uq, we can easily show that Eq. (^E^ can be rewritten 
in a purely thermal form 



Pc 



Tr[e-^^=] ■ 

Making use of Eq. (I2T]) . we find the temperature of the thermal state is 



T 



2arctanh(5/A) 27ra 



(27) 



(28) 



This is exactly the Gibbons- Hawking temperature. Thus, regardless of its initial state, a 
freely falling two-level detector in de Sitter space-time is asymptotically driven to a thermal 
state at the Gibbons-Hawking temperature. Therefore, there must exist a bath of thermal 
radiation in de Sitter space-time. Our open system approach thus reveals that the existence 
of the Gibbons-Hawking effect is simply a manifestation of thermalization phenomena in 
the framework of open system dynamics. At this point, it is worth noting that the Gibbons- 
Hawking temperature of de Sitter space-time has also be derived in other different contexts 
such as the global embedding approach 22|, |23| and the universal Rindler scheme 24 . 



Further aspects of the Gibbons-Hawking effect in terms of the thermalization phenomena 
can be studied by examining the behavior of the finite-time solution (l25ll . For this purpose, 
let us note that 



-2-Ht 



+ 4C2 



|e-4^-Ai + 2e-2(2A+c)r 



A2 cos{Qt) + A; 



sin(f2r) 



(29) 



where 



02 

Ai = [{2A + C)2 + — ] J - 2{2A + C)n + 



A2 
A3 



-2A{A + C)I + {2A + C)U - -U^ , 

2A[— - C{2A + C)]I+ [C{AA + C) - —]n - CV? 



(30) 



Equation (129|) reveals that a freely falling detector in de Sitter space-time is subjected to the 
effects of decoherence and dissipation characterized by the exponentially decaying factors 
involving the real parts of the eigenvalues of 1-L and oscillating terms associated with the 
imaginary part. Therefore, the Gibbons-Hawking effect as a manifestation of thermalization 
phenomena in the framework of open quantum systems actually involves phenomena of 
decoherence and dissipation. This suggests that the vacuum in de Sitter space-time behaves 



like a fluctuating as well as a dissipative medium 



25 



26 



In this regard, our approach to the 
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derivation of the Gibbons-Hawking effect seems to shed new hght on the issue as compared 
to other traditional treatments, and it ties its existence to the effects of decoherence and 
dissipation in open quantum systems . 

In summary, we have analyzed, using the well-known techniques in the study of open 
quantum systems, the time evolution of a freely falling detector in de Sitter space-time in 
weak interaction with fluctuating vacuum conformal scalar fields in the de Sitter-invariant 
vacuum. The detector has been shown to be asymptotically driven to a thermal state at the 
Gibbons- Hawking temperature, regardless of its initial state. Our open system approach 
to the issue therefore shows that the Gibbons-Hawking effect of de Sitter space-time can 
be understood as a manifestation of thermalization phenomena in the framework of open 
quantum systems, which actually involves the effects of decoherence and dissipation. It is 
worthwhile to note that the general techniques developed in the theory of open quantum 
systems may also be applicable to studying other phenomena in curved space-times, such as 
particle creation, and may thus provide new insights in the physical understanding of these 
phenomena. 
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